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Abstract 

An approach to control tethered wings for airborne wind energy is proposed. A fixed length of the lines is considered, 
and the aim of the control system is to obtain figure-eight crosswind trajectories. The proposed technique is based on the 
notion of the wing's "velocity angle" and, in contrast with most existing approaches, it does not require a measurement of 
the wind speed or of the effective wind at the wing's location. Moreover, the proposed approach features few parameters, 
whose effects on the system's behavior are very intuitive, hence simplifying tuning procedures. A simplified model of 
the steering dynamics of the wing is derived from first-principle laws, compared with experimental data and used for the 
control design. The control algorithm is divided into a low-level loop for the velocity angle and a high-level guidance 
strategy to achieve the desired flight patterns. The robustness of the inner loop is verified analytically, and the overall 
control system is tested experimentally on a small-scale prototype, with varying wind conditions and using different 
wings. 

1 Introduction 

Airborne wind energy systems aim at harnessing the wind blowing up to 1000 m above the ground, using tethered wings 
or aircrafts. In recent years, an increasing number of researchers in academia and industry started to investigate this idea 
and to develop concepts of airborne wind energy generators, see e.g. ITT91 ITl [2l [3l [4l [9l ITTl fTOl f8l [23ll as well as lfT3ll 
for an overview. 

In several concepts of airborne wind energy generators that are currently being developed, a tethered flexible wing is 
controlled to fly fast in crosswind conditions, i.e. roughly perpendicular to the wind flow [18 |, and the traction forces 
acting on the lines are converted into electricity using mechanical and electrical equipments installed on the ground 
||2T1 [3l l4l [8l. In particular, on the ground the wing's lines are wound around one or more winches, linked to electric 
generators. Energy is obtained by continuously performing a two-phase cycle, composed by a traction phase, during 
which the lines are unrolled under high traction forces and the generators, driven by the rotation of the winches, produce 
electricity, and by a subsequent passive phase, when the electric generators act as motors, spending a fraction of the 
previously generated energy to recoil the lines. 

The automatic control of the wing is a key aspect of airborne wind energy. In traction phases, the aim is to make the wing 
fly along figure-eight paths, which yield the highest traction forces while preventing line twisting. This control problem 
involves fast, nonlinear, unstable time-varying dynamics subject to hard operational constraints and external disturbances. 
Several contributions by various research groups and companies worldwide have recently appeared in the literature, see 
e.g. H71 |24j [TOj [T31 [TT] [8). Most of the presented approaches ||9] [17] f24] [TO] QUI are based on a nonlinear point-mass 
model of the system, derived on the basis of first-principle laws of mechanics and aerodynamics, and they rely on the use 
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of advanced nonlinear control design techniques. Several proposed techniques employ a reference trajectory derived off- 
line on the basis of the considered model, usually computed in order to be optimal in terms of generated power lfP7] 124] l8l . 
Then, approaches like tracking Model Predictive Control (MPC) ifTTIl or adaptive control-Lyapunov techniques [8] are 
used to track the reference path. The use of MPC with an economic cost function, i.e. without the need to pre-compute 
a reference path, has been also proposed, see e.g. [10], again exploiting a point mass model. While the mentioned 
contributions represent fascinating applications of constrained and optimal nonlinear control methods, their use in a real 
system appears to be not trivial, due to the discrepancies between the employed simplified model and the real dynamics 
of flexible wings, the need to solve complex nonlinear optimization problems in real-time, finally the need to measure the 
wind speed and direction at the wing's flying altitude. 

In contrast with the mentioned works, in a recent contribution, concerned with the control of large kites for seagoing 
vessels IfTTI . a simpler dynamical model has been proposed and used for control design. Such a model has been justified 
by means of measured data, and the designed control system has been tested experimentally, thus showing the practical 
applicability of the approach. It has to be noted that a similar model has been considered also in [8], where it has been 
justified by a priori assumptions. Hence, neither (|8] nor iffD provided an explicit link between the model considered in 
the control design and the wing's characteristics, like area, efficiency and mass. The control approach proposed in IfTTI 
is composed by an inner control loop that computes the input needed to obtain a desired reference heading of the wing, 
and an outer control loop that computes the reference heading according to a bang-bang-like strategy. The inner control 
loop is a quite sophisticated model following approach and it needs the measure of the effective wind speed at the wing's 
altitude, obtained from an onboard anemometer. 

In the described context, we present here new contributions in the field of control of flexible tethered wings for airborne 
wind energy systems with ground-level generators. We focus on the problem of controlling the wing in order to fly along 
figure-eight paths in crosswind conditions. First, we consider a simplified model based on the notion of "velocity angle" 
of the wing, similar in form to the one proposed by [8] and Hill , and we derive an explicit link between the model's 
parameters and the system's characteristics. We show the validity of such a model as compared to experimental data 
collected with a small-scale prototype. This result provides a definitive assessment of the considered control-oriented 
model for tethered wings, thus bridging the gap between theoretical equations and experimental evidence. As a second 
contribution, we present a new control algorithm for tethered wings, based on the derived simplified model. Differently 
from 1 17 24, 8 1, the approach does not employ pre-computed paths based on a mathematical model of the system, thus 
avoiding issues related to model mismatch and to the actual feasibility of the employed reference path for the (uncertain 
and time-varying) dynamics of a real system. Moreover, an estimate or measure of the wind speed at the wing's location, 
as considered e.g. in IfTTI l24l [TOl . is not needed in our approach, nor is a measure of the effective wind speed aligned 
with the wing's longitudinal body axis, used in [11], but only a rough estimate of the wind direction with respect to the 
ground. The structure of the controller is similar to the one proposed in [11], but the controllers employed for the inner 
and outer control loops are different. In particular, the inner controller is a simple static gain, while the outer one is given 
by a switching strategy based on the wing's position. The proposed control approach involves few parameters, that can be 
tuned in an intuitive way. By exploiting the above-mentioned results pertaining to the derivation of the control-oriented 
model, we assess the robustness of the inner control loop analytically, against a wide range of operating conditions in 
terms of wind speed and wing characteristics. We then present the experimental results obtained by testing the approach 
on a small-scale prototype realized at the University of California, Santa Barbara, in different wind conditions and using 
different wings. The paper is organized as follows, section [2] describes the considered layout and the derivation of the 
simplified model for the velocity angle dynamics. The control design is presented in section [3] and experimental results 
are given in section|4] Finally, conclusions are drawn in section[5] 

2 System description and model equations 
2.1 System layout 

We consider a flexible wing, or power kite, connected by three lines to a ground unit (GU). This setup corresponds to a 
prototype built at the University of California, Santa Barbara, shown in Fig. Q] In normal flight conditions, the wing's 
trajectory evolves downwind with respect to the GU. For simplicity, we assume that the nominal wind direction (i.e. 
neglecting turbulence and small, zero-mean deviations) is aligned with the longitudinal symmetry axis of the GU, denoted 
by X. This condition can be achieved by properly orienting the GU, exploiting a measure or estimate of the nominal 
wind direction. With this assumption in mind, our control approach employs the feedback of the wing's position relative 
to the GU to obtain crosswind trajectories, i.e. flying paths that are downwind and symmetric with respect to X axis (i.e. 
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Figure 1: Small-scale prototype built at the University of California, Santa Barbara, to study the control of tethered wings 
for airborne wind energy. 

the wind direction). From our experimental results, it turns out that misalignments of about ±30° between the wind and 
the GU do not change significantly the obtained paths relative to the GU, but indeed give place to lower forces on the 
lines, as expected from the theory (see e.g. [15]). The X axis, together with the Z axis being perpendicular to the ground 
and pointing upwards and with the Y axis to complete a right-handed system, forms the inertial frame G = (X, Y, Z), 
centered at the GU (see Fig. By considering a fixed length of the lines, denoted by r, the wing's trajectory is thus 
confined on a quarter sphere, given by the intersection of a sphere of radius r centered at the GU's location and the planes 
(x, y, z) E M 3 : x > and z > 0. Such a quarter sphere is commonly named "wind window", see Fig. |2](dashed lines). 



4 Z 




Figure 2: Reference system G = (X, Y, Z), wind window (dashed lines), variables 9, <fi, and local north, east and down 
(Ln,Le,Ld) axes. 

The two lateral lines, named steering lines, are linked to the back tips of the wing (see Fig. [TJ and they are used to 
influence its trajectory: a shorter left steering line with respect to the right one impresses a left turn to the wing (i.e. a 
counter-clockwise turn as seen from the GU), and vice-versa. The center line, named power line, splits into two lines 
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connected to the front of the wing (or leading edge) and sustains about 70% of the generated load. The GU is designed 
to support the forces acting on the lines and it is equipped with actuators, able to achieve a desired difference of steering 
lines' length. In the considered prototype, a single motor, together with a linear motion system (visible in the lower-left 
corner of Fig. [TJ, is able to change the difference of length of the steering lines. In this work, we consider a fixed length of 
the lines of r = 30 m and we focus on the problem of designing a controller able to make the wing fly along "figure-eight" 
crosswind paths, which maximize the generated forces. The use of a fixed lines' length does not limit the significance 
of our results, since the problem of crosswind path control can be decoupled from the problem of controlling the lines' 
reeling. Indeed it has been shown that the optimal operation of ground-based airborne wind energy generators is achieved 
with a constant line speed (equal to approximately one third of the wind speed to maximize power production, see e.g. 
lfT31l ). The settings considered here can be seen as a particular case of constant line speed, equal to zero. Other approaches 
in the literature follow a similar strategy of using two separate control systems, one for the wing guidance and one for the 
line unrolling (see e.g. [ 10 8 1 1 1). Rather, the use of such short lines makes the control problem more challenging, since 
the flight paths must be contained in a small area, in order to avoid contact with the ground or aerodynamic stall occurring 
at the border of the wind window. 

The wing is equipped with onboard sensors and a radio transmitter; the receiver and other sensors are installed on the 
GU. The available sensors, together with suitable filtering algorithms, provide accurate real-time estimates of the wing's 
position and velocity vector, to be used for feedback control. For the sake of simplicity, in this paper the feedback variables 
are considered to be measured with no measurement noise; the interested reader is referred to [12] for details on the design 
of suitable filtering algorithms for the described setup. In the next section, we briefly recall a point-mass dynamical model 
of the system and introduce the notion of velocity angle of the wing. The latter represents one of the feedback variables 
in our control approach. 

2.2 Model equations 

We consider the point-mass model already employed in previous works (see e.g. ifTTl [TOl l8l fl4l and references therein) 
and we adapt it to a slightly different reference system, in order to introduce the wing's velocity angle. The latter is one 
of the feedback variables we use and it represents a novelty with respect to the mentioned previous approaches. For the 
sake of completeness, we recall here the main equations of the point-mass model, since they are instrumental to prove a 
theoretical result concerned with the steering dynamics of tethered wings. 

By considering the fixed line length r, the wing's position can be expressed in the inertial frame G by using the spherical 
coordinates 0(t), cf>(t) as (see Fig. [2j: 



where t is the continuous time variable. In (Q]i and throughout the paper, the subscript letter in front of vectors (e.g. cp(t)) 
denotes the reference system considered to express the vector components. 

We define also a non-inertial coordinate system L = (Ln,Le, Ld), centered at the wing's position (also depicted in Fig. 
|2). The Ln axis, or local north, is tangent to the wind window and points towards its zenith. The axis, called local 
down, points to the center of the sphere (i.e. the GU), hence it is perpendicular to the tangent plane to the wind window at 
the wing's location. The Le axis, named local east, forms a right hand system and spans the tangent plane together with 
Lyv- We note that the system L is a function of the wing's position only, and it is different from the local systems used 
in previous works (see e.g. ifTOl and the references therein), due to the different definition of angle 8. A vector in the L 
system can be expressed in the G system by means of the following rotation matrix: 



GP(t) 



rcos (<f){t)) cos (0(t)) 

rsin(0(t))cos (#(*)) 
rsin (9(t)) 



(1) 



R = 




By applying Newton's law of motion to the wing in the reference L we obtain: 



(2a) 



&*) = L t L j\! + 2tan(0(i))0(t)0(t), 




(2b) 
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where m is the mass of the wing. In (f2ab-d2bl> and throughout the paper we denote unit vectors by "e" followed by a 
subscript indicating the related axis, e.g. eL N (t) denotes the unit vector of the ijy axis. The force LF(t) consists of 
contributions from the gravity force hF g (t) and aerodynamic force lF z (£). Vector LF g (t) can be computed as: 



LF g (t) = 



— mg cos (0(t)) 


mg sin (6(t)) 



(3) 



where g is the gravity acceleration. The aerodynamic force is given by the contributions of the lift and drag generated by 
the wing and of the drag induced by the cable. These forces depend on the the effective wind, W e (t), computed as: 



W e (t)=W(t)-v(t), 



(4) 



where W(t) is the wind relative to the ground and v(t) = -^p(t) is the wing velocity vector, which can be expressed in 
the L frame as 



At) 



r6{t) 
rcos ((9(i))0(t) 




Then, the aerodynamic force LF a (t) can be computed as (see e.g. ifTOl ): 

FS) =\pC L {t)A\W e {t)\ 2 z w (t) + 



-pC D {t)A\W e {t)\ 2 x w {t)+ 
^ P C D jA[ cos{Aa(t))\W e (t)\ 2 x w (t) 

±pC L {t)A\W e {t)\ 2 z w {t) + 

1 (n m i C pjA cos (Aa(t)) \ 2 ^ 
-p \C D {t) + — I A\W e (t)\ x w (t). 



(5) 

(6a) 
(6b) 
(6c) 

(6d) 



C D ,e,(t) 

In (0, the contributions d6at-(l6bl are, respectively, the lift and drag forces generated by the wing, while (l6ct is the drag 
induced by the lines. C^it) and Coit) are the aerodynamic lift and drag coefficients of the wing, Cd,i is the drag co- 
efficient of the lines, A is the reference area of the wing, Ai is the reference area of the lines, p is the air density, and 
x w (t) and z w (t) are the directions of the drag and lift forces, respectively. The parameter CD,eq{t) is called the equivalent 
aerodynamic drag coefficient, since it accounts for the drag of both the wing and the lines. We note that the aerodynamic 
coefficients are considered as time-varying parameters here, since they depend on the wing's angle of attack, which in 
turns changes in time as a function of the flight conditions. The variable Aa(t) is the angle between the effective wind 
vector W e (t) and the tangent plane to the wind window at the wing's location. The vectors x w (t) and z w (t), defining the 
directions of the lift and drag forces, depend on the direction of the effective wind and on the roll angle ijj(t) of the wing. 
In particular, x w (t) points in the direction of the effective wind W e (t), while z w (t) is perpendicular to x w {t) and to a 
further vector, denoted by et(t), which points from the right tip of the wing to the left one, as seen from the GU (see e.g. 
IfTUl for a formal definition). 



Vectors x w (t) and z w {t) can be expressed in the L frame as: 

LX w (t) = 



-COS (£(*)) 

-Bin (£(«)) 




-sin (*(«)) 
cos (£(*)) 
-1 



LZ w (t) = 



(cos (Aa(t)p 

[sin (Aa(t)), 

-con (£(*)) -sin(£(i)) 0\ 
-sin(CW) cos (£(<)) 
-1/ 



(7a) 



(7b) 



' -cos(?/>(t))cos(7;(t))sin(Aa(t)) N 

cos (i/j(t)) sin sin (Aa(t)) + sin (ip(t)) cos (Aa(t)) 

i cos (ip(t)) cos (rj(i)) cos (Aa(t)) , 
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In (0, rj(t) is given by (see e.g. 171): 



rj(t) = arcsin (tan {Aa(t)) tan (V> (*))), (8) 
and is a function of the steering input, 5(t): 

ip(t) = arcsin (9) 

where d s is the wing span. Finally, is the heading angle of the wing, and it is computed as the angle between the local 
north ijy and the effective wind W e (t) projected on the (Ln,Le) plane: 

^) = arctanf^^V (10) 
\W e {t)-e LN (t)) 

In dTOb . the four-quadrant version of the arc tangent function shall be used, such that G [— tt, it]. The assumption 
underlying equation dTOb for the computation of the wing's heading is that the wing's longitudinal body axis is always 
contained in the plane spanned by vectors W e (t) and p(t). 

Equations (f2ab-(fT0b give an analytic expression for the point-mass model of the wing, with four states (6(t), 4>(t), 9(t), <j>(t)), 
one manipulated input (S(t)) and three exogenous inputs (the components of vector W(t)). Such a model has been widely 
used in the literature on control design for airborne wind energy applications, see e.g. 10 [17] [10] [8] . Several existing 
approaches are based on nonlinear model predictive control techniques |[T7l[T0l , leading to quite complex multivariable 
controllers that rely on the feedback of the four states. In the section l2~4l we show how a simpler model, which we will 
use for our control design, can be derived from the dynamics (T2ab-(fT0b. The variable involved in such control-oriented 
model is the velocity angle "f(t) of the wing, defined as: 

( v{t) ■ e L Jt)\ f cos (6(t))<p(t)\ 

-lit) = arctan _) { ) { = arctan v , v nvy ' . (11) 

y> \v{t) ■ ?L N {t) J \ 0(f) J 

The angle 7(f) is thus the angle between the local north eL N (t) and the wing's velocity vector v(t). This variable is 
particularly suited for feedback control, since it describes the flight conditions of the wing with just one scalar: as an 
example, if 7 = the wing is moving upwards towards the zenith of sphere, if 7 = tt/2 the wing is moving parallel to the 
ground towards the local east, finally if 7 = it the wing is flying towards the ground. Hence, the time derivative 7 defines 
how fast the wing is being steered while flying in the wind window. Similarly to (fT0T >, also in ( fTTT i the four-quadrant 
version of the arc tangent function shall be used, such that 7(f) E [— tt, it]. 



2.3 Input model 

The prototype used for our test flights features two attachment points for the steering lines on the GU, left and right. These 
attachment points are separated by a distance d (see Fig. [1] in the lower part, where the attachment points with swaying 
pulleys are visible). When the wing's lines are not aligned with the X axis, this distance induces an equivalent steering 
deviation. We call such deviation "geometric input", S g , since its value depends on the geometry of the attachment points 
and on the (9, (f>) position of the wing in the wind window. Hence, the overall steering input acting on the wing is: 

5{t) = 6 a (t) + <y g (t), (12) 

where S u (t) is the input issued by the control system, i.e. the difference of length of the steering lines (right minus left) 
obtained by changing the position of the linear motion system on the GU. In this section, we derive an expression to 
compute the geometric input as a function of 9[t), <j>{£) and d. Consider the GU and the wing as seen from above, such 
that the trace of the wind window at the wing's height is a semicircle of radius r cos (6) and, for fixed 9, the wing position 
is univocally identified by <f>. Let us consider first a situation in which both wing tips lie on the tangent plane to the wind 
window at the current wing position (6, (j>). Since the steering of the wing is essentially induced by a roll motion, by 
which the wing tips are moved away from the tangent plane, we call this orientation "neutral configuration" (see Fig. [3] 
gray drawing). Now, assuming that S a = 0, the left and right lines have the same length, so that the wing tips are forced 
to leave the neutral configuration (see Fig. [3] black drawing), in the same way as a steering input were acting on the wing. 
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In other words, in order for the wing to be in a neutral configuration in the presence of a value of <fi ^ 0, the two steering 
lines should have a difference of length equal to d sin (0) cos (9). Hence, we can compute the geometric input as: 



*g(*) 



d sin (<£(£)) cos (6(t)). 



(13) 



We note that the geometric input is always bounded by d, moreover with an increasing 9 position of the wing its value 
decreases and becomes eventually zero if the wing is at the zenith position of the wind window. The minus sign in ( fT3l 



Figure 3: Sketch of the wing in a generic (9, 4>) position as seen from above, and related geometric input S g = 

—d sin {(f)) cos (9). 

derives from the fact that the geometric input imposes a counter-clockwise turn (as seen from the GU) for <f> < (i.e. 
7 > 0) and a clockwise turn for <j> > (i.e. 7 < 0). 

2.4 Control-oriented model for tethered wings 

In ifTTI . a simple model was presented and used for the control design, where the time derivative of the heading angle 
of the wing is given as a function of the control variable. Here, we consider a similar simplified model, where we use 
the velocity angle instead of the heading angle. This choice is supported by the assumption of small sideslip angle, as 
formally stated below: 

Assumption 1 The difference between the velocity angle 7(4) and the heading angle is negligible, i.e. the effective 
wind projected onto the tangent plane to the wind window at the wing's location is equal to the wing's velocity v(t). 
Moreover, all the forces in the direction of vector v(t) are negligible as compared to lift and drag. M 

Assumption [JJ is common in the analysis of airborne wind energy generators (see, e.g., lfl"8] [151 ) and it is reasonable 
whenever the wing is flying downwind roughly perpendicularly to the wind flow. As a further justification for this 
assumption, we show in Fig. [4] the good matching between the values of j(t) and measured during figure-eight 
paths carried out with our small-scale prototype and a 9-m 2 power kite. With Assumption[TJin mind, the model proposed 
in IfTTI can be written in the following form: 



In particular, in IfTTI the term K (t) is given by a constant gain, qk, multiplied by the magnitude of the effective wind speed 
aligned with the wing's heading. Such a model was justified mainly through experimental results in [ 1 1 1, where the gain 
gx was derived empirically. In J8], a model of the form ([14) is also used, and justified by some assumptions. Equation 
(fl4l > provides indeed a model well-suited for control design, with the right balance between accuracy and simplicity, 
however neither IfTTI nor JU derived an explicit relationship between the model's parameters and the main characteristics 
of the system, such as wing size, mass, or efficiency. The theoretical result we are presenting next is aimed to fill this gap, 
by linking equation (fT4l > to the first-principle model recalled in section [2721 Moreover, we present experimental data that 
confirm the validity of our result. We consider the following assumption on the roll angle i/j(t): 




y(t)~K(t)S(t)+T(t). 



(14) 
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Figure 4: Experimental results. Angles 7 (solid line) and £ (dotted) measured during controlled flights in crosswind 
conditions. Wing size: 9-m 2 . 



Assumption 2 The angle is sufficiently small to linearize its trigonometric functions. M 

Assumption|2]is reasonable in the considered context, since for example the prototype at the UC Santa Barbara operates 
with tp ~ ±7.5° (for a 12-m 2 wing with wingspan d s — 3.1m) up to ip ~ ±12.8° (for a 6-m 2 wing with wingspan 
d s = 1.8 m). 

We can now state our theoretical result. 

Proposition 1 Let assumptions\T§2\hold. Then, equation (1141) holds with: 



pC L (t)A / , 1 



2 



t = sa >w))^ Mt )) + sta k 

K*)l 

where \v(t)\ is the magnitude of the wing's velocity, and E eq (t) = C / 'C 'd ,eq(t). 
Proof 1 See the Appendix. 

Proposition Q] provides an explicit link between the main lumped parameters of the wing, like area, mass and lift coef- 
ficient, and the gain and external disturbance of (fl4l . It is worth elaborating more on this result and its implications. 
According to equation (fl4l . there is basically an integrator between the control input S(t) and the velocity angle 7(f), 
with a time-varying gain. In the first term of (I15at . we can see that such a gain increases as the wing's speed does; thus, 
a larger speed provides higher control authority but it can also bring forth stability issues, if the control system is not 

properly designed. For a given wing flying in crosswind conditions, the term 9 i^ff - ( 1 + (t) ) c ' oes not cnan g e 
significantly during operation, hence supporting the model proposed by [11], where, as mentioned above, a constant gain 
qk , multiplied by the effective wind speed, is used to compute K. Equation d!5ab also implies that a larger area-to-mass 
ratio gives in general a higher gain K, and that the steering behaviors of wings with similar design (i.e. similar aero- 
dynamic coefficients) but different sizes are expected to be similar, provided that the area-to-mass ratio, A/m, does not 
change much when the size is scaled up. However, the latter consideration holds true if just the wing itself is considered, 
since for example the use of larger wings would require lines with larger diameter, and the consequent added mass and 
induced drag would generally reduce both the equivalent efficiency and the area-to-mass ratio. 

Furthermore, since the wing's speed is roughly proportional to the effective wind speed projected along the lines' direction 
(see e.g. lfT31l ). the gain K can be also re-written as (assuming the wind is aligned with the AT-axis, as we remarked in 
section l2~TT i: 

k { t) = pCL(t)AE^t)cc B m))ccs(4>(t ) ) a 1 y m 

2md s V Elqit)) 
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(c) 



(d) 





Figure 5: Experimental results, (a)-(c): Comparison between the measured values of 7(t)/ I 1 



El q {t) 



as a 



pC L A 
2md s 



as 



function of the steering input <5(t) (gray and black dots) and the theoretical linear relationship given by the gain 
per equation d!5ab (solid line). The gray dots represent experimental data collected in the whole range of 6, tj> spanned by 
the wing during operation, while the black dots represent values collected when \<p\ < 5°, i.e. in crosswind conditions, (a) 
Airush One® 6 kite, (b) Airush One® 9 kite, (a) Airush One® 12 kite. The lumped parameters for the kites are reported 
in TableQ] (d): Comparison between the value of 7 obtained in experimental tests (solid) and the one estimated using the 
simplified model and the result of Proposition Q] (dashed). Wing size: 6-m 2 . 



According to (TToT i. the value of K is highest when the wing is flying crosswind, and it decreases as the wing approaches the 
borders of the wind window (i.e. 9 ~ and/or <f> ~ ±§). Therefore, maneuvering the wing in these conditions requires 
larger control inputs and in some situations it might be not possible, leading to a complete loss of controllability. Finally, 
the wing's efficiency and lift coefficient also play an important role: in particular, for fixed drag coefficient the gain is 
expected to grow quadratically with the lift coefficient, as it can be noted from (TToT i by considering that E eq — Cl / Co,eq- 
The term T(t) basically accounts for the steering effect that the gravity and apparent forces have on the wing. It can be 
noted that the influence of gravity gets smaller with larger wing speed, hence it can be easily dominated by the control 
action when the wing's lines are aligned with the wind, while it becomes more and more important as the wing moves to the 
side of the wind window. The result of Proposition[T]is also well-confirmed by the experimental tests we carried out with 
our small-scale prototype. In order to show the matching between our result and the experimental evidence, we compare 

the measured values of 7, normalized by ^1 + E -i~ ^ J W(t)\> an d the steering input 5(t): according to equation jl5ab . 

these two quantities are proportional through the gain given by p ^ d A ■ We show the matching between the experimental 
data and such a relationship in Fig. |3a), (b) and (c), for a 6-m 2 , a 9-m 2 and a 12-m 2 wing, respectively. The lumped 
parameters of the employed wings are resumed in Table Q] In Fig. |3a)-(c), the gray dots represent experimental data 



9 



collected in the whole range of 6, 4> spanned by the wing during operation, while the black dots represent values collected 
when < 5°, i.e. in crosswind conditions. It can be noted that the linear relationship computed by using the lumped 
parameters as given by the model (fT5t matches quite well with the experimental data, not only in crosswind conditions, 
where the underlying assumptions are valid, but also with larger values of tj>, in the range ±35°. Finally, in Fig. |5}(d) 
we show an example of the matching between the time course of 7 during figure-eight paths with the 6-m 2 wing and 
the estimate given by the model (TT4-b - (TT3T> . Overall, the presented results provide a solid bridge between the existing 
dynamical models of the steering behavior of tethered wings and experimental evidence. In the next section, we introduce 
a new control approach for crosswind flight of tethered wings, which is based on the presented simplified model. 

Table 1: Lumped parameters of the wings employed in the experimental activities. 



Airush One® 12 


Area 


A 


12 m 2 


Mass 


m 


2.9 kg 


Wingspan 


d s 


3.1m 


Lift coefficient (average) 


C L 


0.85 


Equivalent aerodynamic efficiency (average) 


E eq 


5.3 


Airush One® 9 


Area 


A 


9 m 2 


Mass 


m 


2.45 kg 


Wingspan 


d s 


2.7 m 


Lift coefficient (average) 


C L 


0.8 


Equivalent aerodynamic efficiency (average) 


E eq 


5.6 


Airush One® 6 


Area 


A 


6 m 2 


Mass 


m 


1.7 kg 


Wingspan 


d s 


1.8m 


Lift coefficient (average) 


C L 


0.6 


Equivalent aerodynamic efficiency (average) 


E eq 


5.1 



3 Control design 

We propose a control scheme consisting of three nested loops, shown in Fig. [6] The outer control loop employs the current 
wing position, in terms of 8, <fi angles, to compute a reference velocity angle, 7 re f, for the middle control loop. The latter 
employs 7 as feedback variable and it has the objective of tracking 7 re f, by setting a suitable position reference, <5 mjre f, 
for the actuator installed on the GU. The innermost control loop then employs a feedback of the motor position S m to 
command the motor's current i m , in order to track the desired position S mjs f . This control structure allows to separate the 
nonlinear part of the controller, which is all kept at the outermost level, from the linear one, hence obtaining two simple 
controllers for the middle and inner loops, for which we can carry out a theoretical robustness analysis. The control 
algorithm for the outer loop is given by a quite simple guidance strategy, hence making the whole control system very 
suitable for implementation and experimental testing. In the next sections, we describe each control loop in details. We 
first adopt a general notation for the involved design parameters and we then provide in section |4]the specific numerical 
values employed in our tests. 
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Figure 6: Overview of the proposed control system. 



3.1 Position control loop 



The innermost loop consists of a standard position control system with an electrical DC brushed motor and a linear motion 
system, based on a lead screw mechanism. By neglecting high-order effects, in the absence of external disturbances the 
dynamics of the actuator can be modeled as: 



Lit) 



(t)+K m i m (t), 



(17) 



where S m (t) is the actuator's position in m, K m , uj m are parameters depending on the actuator's characteristics and i m (t) 
is the commanded current in A. Moreover, the prototype is equipped with a series of pulleys such that, for a given value 
of S m , the corresponding difference of length of the steering lines is equal to 



5 u (t) = K s 6 m (t), 



(18) 



with Kg being a constant gain. The value of 6 m is measured with high accuracy by an optical rotary incremental encoder, 
suitably scaled to obtain the linear position of the actuator from the motor's angular position. We use standard loop- 
shaping control design techniques (see e.g. Il20l ) to design a cascade feedback controller C for this control level. The 
control scheme is shown in Fig. [7] We note that the actuator has current limits of ±i m , included as a saturation in Fig. [7] 



4vef + 



c 




I 




Actuator 













Figure 7: Scheme of the actuator position control loop. 



Such a saturation does not give rise to integrator windup problems, since the employed controller C does not have integral 
action. Yet, the presence of an integrator in the plant (IPTl still yields a closed-loop transfer function with unitary gain 
between the reference and the output. In particular, the closed loop system for the inner loop results to be of the form: 

<L(*) + 2( cl UjJ m (t) + L0l5 m {t) = W c Vm,ref(t), (19) 

where <5 m , re f (t) is the reference position provided by the controller for the wing's velocity angle, which is described next. 

3.2 Velocity angle control loop 

The control design for the middle loop exploits the control- oriented model (fl5l > derived in section l2~4l The control input is 
the actuator position <5 m ,ref ) (i-e- the reference for the innermost control loop described in section lXTT ). and the feedback 
variable is the wing's velocity angle ^(t). The controller is given by a simple proportional law : 

<WW = K c (7ref (t) - 7tt)) , (20) 
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where 7 re f (t) is the target velocity angle provided by the outer control loop, and K c is a scalar gain to be chosen by the 
designer. The value of K c is the only design parameter for the velocity angle controller, and it can be tuned at first by using 
the equations dT4b-(fT~5b with some nominal system parameters, and then via experiments. Moreover we note that, with the 
proposed approach, K c can be tuned in order to robustly stabilize the control loop. In particular, for a fixed value of K c , 
by exploiting the model $15[ we can carry out a robustness analysis of the control system comprising the velocity angle 
loop and the actuator position loop. Such analysis allows us to discern for what range of wind speeds and wing's parame- 
ters the considered static gain K c is able to stabilize the velocity angle control system. In order to do so, we re-write the 
dynamical system given by equations (ITSb-dTSb and (TT9b-(f20T> in terms of velocity angle tracking error, A 7 (t) == 7 re f— 7(£): 



(21) 



" A T (t) ■ 




-KKs 




■ A 7 (t) ■ 




L(t) 




1 




<5m(*) 


+ w(t) 


L(t) 













A d {K) 

In (fJTJ, the term K corresponds to the gain in (fT5t and depends on the system's parameters as well as the wind and 
flight conditions. The term w(t) accounts for the effects of gravity and apparent forces of (TT31 ), as well as the forces 
exerted by the lines on the innermost control system, and it can be seen as a bounded disturbance acting on the sys- 
tem. System $2l[ has time-varying, uncertain linear dynamics characterized by the matrix A C \(K). In fact, the scalar 
K(\v(t)\, A, E eq (t),CL(t),m,d s ) is a function of several parameters, like the wing's speed magnitude and the wing's 
efficiency and lift coefficient, that are not precisely known and vary over time, since for example the aerodynamic coef- 
ficients depend on the wing's angle of attack, which changes in time according to the flight conditions. However, upper 
and lower bounds for all of the involved parameters can be easily derived on the basis of the available knowledge on the 
system, and these bounds can be employed to compute limits K\ i = 1, 2, such that K G [K 1 , K 2 ]. Then, any possible 
matrix A C \(K) results to be contained in the convex hull defined by the vertices A^K 1 ), A C \(K 2 ). After computing these 
vertices, existing results based on quadratic stability and linear matrix inequalities (LMI) can be used to assess the robust 
stability of system (f2Tb . In particular, system (fJT) results to be robustly stable if there exists a positive definite matrix 



P = P 2 € 



5x3 



such that (see e.g. |6|): 



Al{K l )P + PA c i(K l ) -< 0, i = 1, 2, 



(22) 



where T stands for the matrix transpose operation. Condition (l22t can be easily checked by using an LMI solver. In 
section 0] we show that indeed a unique value of K c guarantees robust stability of the velocity angle control loop for a 
wide range of operating conditions. Since the wing's velocity can be measured quite accurately, one could also adopt a 
gain scheduling approach and make the gain in ( f20b depend on \v(t)\, in order to improve the performance. However, in 
our experimental tests this was not needed, as the control systems achieved satisfactory performance for all the conditions 
of wind speed (and, consequently, of wing speeds) that we experienced, as well as with all the three different wings we 
tested. 

Finally we remark that, due to the physical limitations of the actuator, the position S m (t) is constrained in the range ±S m , 
where S m is a positive scalar. Similarly, we saturate the reference position 5 m<le f (t) in the same range. However, we note 
that these saturations are never active during operation (we show this fact with measured data in the next section), hence 
the stability analysis reported above is still valid even if it does not account explicitly for the input limits. Fig. |8]shows a 
scheme of the velocity angle controller 



K. 



m,ref 



Figure 8: Scheme of the velocity angle controller. 



3.3 Outer control loop 

The outer control loop (shown in Fig. [9} is responsible for providing the velocity angle control loop with a reference 
heading j re f{t). The goal is to have a control algorithm able to make the wing fly along figure-eight paths, with few 
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Figure 9: Scheme of the guidance strategy. 



tuning parameters. Hence, we want to avoid the need to pre-compute a whole trajectory to be used as reference, as it has 
been done in previous works |[T7l l8l. In fact, pre-computed reference paths are based on some mathematical model of the 
system, with the consequent unavoidable issues of model mismatch and approximation that might give rise to problems 
related to stability and attractiveness of the chosen trajectory. In ITTT1 . a bang-bang like strategy to set the reference 
heading of the wing, avoiding the use of pre-computed reference flying paths, has been described. Here, we propose a 
similar approach, whose advantage is to provide a simple and explicit link between the tuning parameters and the position 
of the resulting paths in the wind window. 

In our approach we define two fixed reference points in the (<j), 9) plane, denoted by P_ = (</>_ , 0_) and P + = ((/)+, 9 + ), 
with 0_ < <j) + (see Fig. ITOb . The controller computes a new value of the reference velocity angle at discrete time instants. 





— r — pw7~Wwing\ 
v- \ \ 
<L \ \ 

\ \ 


1 Active Target Point ^ 


i 

i 





Figure 10: Sketch of the control strategy for the outermost loop. Wind window projected on the (Y, Z) plane (black solid 
line), target points P_ , P + ('o'), traces of points with constant 9 and (f> (dashed lines), example of wing's path (gray solid 
line), and example of how the reference 7 re f is computed for a given wing position ('•'). 

At each time step k G Z, one of the two reference points is set as the active target P a (k) = (<f> a (k), 9- d (k)), according 
to a switching strategy that we describe next. Then, a first target velocity angle j re f(k) is computed on the basis of the 
measured values of 9(k) and </>(fc) as follows: 

i.e. in order to make the wing's velocity vector point towards the active target (compare the definition of j(t) in (fTTli). 
Finally, the actual reference 7 re f is computed as the output of a 2 nd -order Butterworth filter F^, whose input is 7, e f, with 
cutoff frequency cj 7 , to be tuned by the control designer, in order to provide a smooth reference as input to the velocity 
angle control loop. The target points are switched according to the following strategy: 

M(j>{k) < 4>- then P a (fc) = P+ 

If <j>{k) > <p+ then P a (fc) = P_ (24) 
ElseP a (fc) = P a (fc - 1). 

Thus, the target point is switched when the measured value of 4> is outside the interval [</>-, 4>+]- By how we defined the 
velocity angle 7, after the target point has been switched the wing will start turning, under the action of the inner control 
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loops, following "up-loops", i.e. pointing first towards the zenith of the wind window and then to the other target point. 
Equations (I23ll-(l24li provide the core of the controller we propose for setting the reference velocity angle. We note that 
also this control loop involves few parameters, i.e. the target points P_, P + , which can be intuitively tuned by using 
simplified equations of crosswind motion like those presented in [ 15 1, and the cutoff frequency w 7 , which can be tuned in 
order to have sharper (for higher w 7 ) or wider (for smaller u> 7 ) turns of the wing (we provide an example obtained from 
our experiments in the next section). Finally, we note that the proposed control strategy does not rely on any measurement 
or estimate of the wind speed, however the interval [</>_ , <p + ] shall be centered around the wind direction, in order to make 
the wing fly in crosswind conditions. A rough measure of the wind direction can be obtained by means of standard wind 
vanes, moreover an estimate can be obtained from other measured quantities, like the wing's speed and line forces. 

4 Experimental results 

We implemented the controller described in section [3] on a real-time machine made by SpeedGoat® and programmed 
with the xPC Target® toolbox for MatLab®. We employed three different wings in our tests: a 6 m 2 , a 9 m 2 and a 12 m 2 
Airush® One power kites. The lumped parameters of the wings are reported in Table [TJ The sampling frequencies we 
used for the control loops are 100 Hz for the innermost controller and 50 Hz for the middle and outermost ones. A movie 
of the experimental tests is available online [5|. Table [2] shows the main system's and control parameters, in addition to 
those reported in Table [TJ In Table [2] the values of the target points P + , P_ and of the cutoff frequency o> 7 are indicated 
as "variable" because we tested different values, in order to assess the influence of these tuning parameters on the obtained 
performance. The first set of results that we show (Figs. n~TTTT~4-b is related to the 9-m 2 wing; the target points' coordinates 
are 9+ = 0_ = 0.35 rad, 0_ = — 0.2rad, 0+ = 0.2 rad and the value of w 7 is 0.25 Hz. 

Fig. Qjjshows an example of the courses of <5 mire f(i) and S m (t) obtained during the experiments with the employed con- 

0.05 1 1 1 1 1 1 1 



0.04 1 




Time (s) 

Figure 11: Experimental results. Example of courses of target position <5 m>re f(t) (dashed line) and actual position 5 m (t) 
(solid line) of the actuator obtained with the designed innermost control loop, during automatic flight tests of a 9-m 2 wing. 
The numbered circles corresponds to the conditions highlighted in Fig. Q~2] Employed wing: 9 m 2 . Guidance parameters: 
6+ = 9- = 0.35 rad, </>_ = -0.2 rad, 0+ = 0.2 rad; w 7 = 0.25 Hz. 

trailer for the innermost control loop. Note that, for the sake of simplicity, we did not include the presence of disturbances 
in the description of this control loop in section [XT] Indeed the forces applied by the wing's lines on the actuator can be 
modeled as an additive disturbance at this level. The effect of such disturbance can be clearly seen between points 8-1 and 
around point 5 in Fig. [TTJ where there is some error between the reference and actual position. By how the machine has 
been designed, the force exerted by each steering line on the actuator is equal to twice the force on the line. This amounts 
to approximately 800 N at the mentioned points in the figure. The tracking error induced by such disturbance does not 
influence much the performance of the overall system. If needed, the position controller can be tuned to achieve better 
tracking performance, at the cost of higher energy consumption. The flown path corresponding to Fig. [TTJcan be seen in 
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Table 2: System and control parameters. 



System's parameters 


Actuator gain 


K m 


0.73m/(sA) 


Actuator pole 




1.9rad/s 


Position limits 


s m 


0.35 m 


Current limits 




1 A A 
10 A 


Distance between 

steering lines' attachment points 


d 


0.5 m 


Tether length 


r 


30 m 


Gain between motor position 
and line lengths' difference 


5u 


4 


Air density 


P 


1.2kg/m 3 


Position control loop 


Sampling frequency 




100 Hz 


Damping 


Ccl 


0.7 


Natural frequency 




78 rad/s 


Velocity angle control loop 


Sampling frequency 




50 Hz 


Feedback gain 


K c 


0.046 m/rad 


Guidance strategy 


Sampling frequency 




50 Hz 


Target points 


P+,P- 


variable 


Butterworth filter 
cutoff frequency 


w 7 


variable 



Fig. [12] as well as the line forces and the velocity angle of the wing. In order to evaluate the robustness of the velocity 
angle control loop, we checked the condition d22l for the following ranges of the involved parameters: \v\ £ [2, 80] m/s, 
E eq e [2, 8], C L £ [0.4, 1], A e [6, 12] m 2 , d s 6 [1.8, 3.1] m, m £ [1.7, 3] kg. The considered ranges of wing speed and 
equivalent efficiency correspond to wind speeds in the interval [1, 10] m/s, which covers a wide range of normal operating 
conditions for airborne wind energy generators. With the considered intervals, the control systems results to be stable with 
the chosen value of the gain K c = 0.046 m/rad, thus indicating a good robustness of the approach. Indeed in our tests the 
same value of K c was used for all wind conditions and all three wings, with good results. We remark that, as anticipated 
in section l3~2l the saturation 8 m = 0.35 m on the motor position was never active since, in practice, the commanded 
position is much lower, of the order of 0.1 m (see Fig. [TTV The wing can be effectively steered with such small values 
of control input thanks to the presence of the geometric input S g , whose effect is to contribute to steer the wing towards 
the center of the wind-window, hence facilitating the desired up-loops figure eights. In other words, the geometric input 
already gives place to a "self-steering" behavior, and the aim of the controller is to apply relatively slight corrections to 
prevent instability and divergence of the flown paths, especially in the middle of the wind window, where the geometric 
input is small and the wing, without feedback control, would head straight towards the ground. This can be clearly seen by 
comparing the numbered points in Figs. QT]and[T2](top). On the bottom of Fig. [12] a typical example of the courses of the 
reference 7 re f(i) and of the actual velocity angle ^(t) obtained during the experiments is shown: it can be noted that the 
middle control loop achieves good performance in tracking the desired velocity angle. The proposed control approach has 
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Figure 12: Experimental results. Single figure-eight path obtained during automatic test flights with about 2.4 m/s wind 
speed. From top to bottom: flying path in (cj>, 8) coordinates, course of the total force acting on the lines (solid line) and 
of the forces acting on the left (dotted), right (dash-dot) and center (dashed) lines, course of the velocity angle 7 (solid 
line) and reference velocity angle 7 re f (dashed). Employed wing: 9 m 2 . Guidance parameters: 8 + = 8- = 0.35 rad, 
<t>- = -0.2 rad, <j)+ = 0.2 rad; u 1 = 0.25 Hz. 



been successfully tested under various conditions. The controller was able to deal with varying wind speed between about 
2 m/s up to 6 m/s, achieving similar, consistent flight paths also in the presence of gusts. Wind speeds lower than 2 m/s 
would not allow the employed wings to fly without stalling, and we avoided to test with wind speeds higher than 6 m/s not 
to stress too much the wing itself, the lines and the mechanical frame, pulleys and other components of the prototype (we 
recall that the involved forces increase linearly with the square of the wind speed). Misalignments of the GU with respect 
to the wind direction up to 30° did not pose a problem for the overall control strategy. However, if the target points (and 
thus the figure-eight paths) are not centered with respect to the wind, the flight trajectory becomes slightly asymmetric 
in terms of altitude. In Fig. Q~3] this phenomenon can be seen in the (<fi, 8) plane. The forces are larger on one side of 
the loop and the wing tends to gain more altitude during these turns. Fig. [14] shows the same paths as Fig. Q~3] but in 
the (X, Y, Z) frame. Fig. [l5]shows the results obtained with different wind speeds, again with the 9-m 2 wing and target 
points set to 8 + = 8- = 0.8, rad, 0_ = —0.4 rad, cj>+ = 0.4 rad. The results of Fig. [I5]highlight the fact that, since the 
proposed controller does not aim to track a given, specific reference path, the resulting figure-eight trajectories change 
with different wind conditions. However, the results also indicate that such changes are not dramatic, moreover the flown 
paths can be intuitively adjusted by changing the position of the target points P_ , P + employed in our control approach. 
As mentioned in section 13.31 the cutoff frequency w 7 can be used to influence the course of the reference velocity an- 
gle. The higher this frequency, the faster the transient behavior of 7 re f(t) when the target point is switched, with con- 
sequent sharper turns and smaller flown paths: an example can be seen in Fig. [16] where the target points are set to 
8 + = 8- = 0.55, 4>- = —0.2, 4> + = 0.2 and two different values of uj 1 are considered. Figs. [P71 and 1X81 show an analysis 
of typical figure-eight paths obtained with the 12-m 2 and with the 6-m 2 wings, respectively: while the qualitative results 
are similar to those obtained with the 9-m 2 wing, these Figures show the quantitative differences in terms of generated 
forces and behavior of the velocity angle. Finally, we present a comparison of the forces generated during extensive ex- 
perimental tests with the different wings, and the corresponding theoretical values obtained e.g. from the results of [ 15 1. 
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Figure 13: Experimental results. Ten consecutive figure-eight paths in the (</>, 8) plane, and corresponding linear dimen- 
sions with the employed line length of 30 m. Employed wing: 9 m 2 . Guidance parameters: 6+ = 6*_ = 0.35 rad, 
6- = -0.2 rad, 4> + = 0.2 rad; w 7 = 0.25 Hz. 



In particular, the generated force is expected to be affine in the quantity E eq yl + -g^-j ( cos 6 cos(<j>)\W\j , with the 

gain being equal to pC J; — , The obtained results, shown in Figs. [T9l a)-(c), indicate that the qualitative behavior is indeed 
consistent with the theory, albeit with some variability. The latter is mainly due to the uncertainty in the measurement of 
the wind speed at the wing's height and in the estimate of the lift coefficient and equivalent efficiency. 



5 Conclusions 

A simplified model for the steering behavior of tethered wings in crosswind flight has been derived, and its validity has 
been assessed with experimental data, collected with three different wings. Then, an approach to design a feedback 
controller for tethered wings has been proposed, with the aim to obtain figure-eight crosswind flying paths, to be used in 
airborne wind energy generators. The controller features three hierarchical levels. Differently from existing approaches 
in the literature, neither the measurement of the wind speed at the wing's altitude nor that of the effective wind speed 
are required. The control system involves few parameters, that can be easily tuned. Moreover, a robustness analysis 
of the inner control loop has been carried out. The effectiveness of the approach has been shown through extensive 
experimental results obtained with a small scale prototype and different wings. Finally, a comparison between the forces 
measured during the tests and the ones predicted by the existing theoretical results has been carried out, showing a general 
qualitative consistency between theory and experiments. 

The natural development of this research will be the execution of full generating cycles with the considered concept of 
airborne wind generator, and the comparison of the obtained results with the existing theoretical and numerical analyses. 
In order to achieve this goal, an additional motor/generator has to be used, and the related control algorithms have to 
be designed. These tasks could not be pursued with the considered prototype, which does not have energy generation 
capability, and they are subject of current research. 

Proof of Proposition^ For the sake of simplicity of notation, in the following we drop the dependance of time-varying variables 
from t. We start by deriving an analytic expression for the velocity angle rate. For this we take the derivative of ill 1 1 to get 

cos (61) 6<j> - sin (61) 4>Q 2 - cos (6>) <j>6 
1 ~ cos (6») 2 4> 2 + e 2 

The accelerations 9 and (f> in \25\ can be expressed as functions of the forces acting on the wing by using the model equations ll2al>-( ll01 >. 
By Assumption Q] neglecting all forces except for the aerodynamic ones and considering the balance of the lift and drag forces in the 
direction of the wing velocity v we have (see e.g. [ 15 1): 

sin (An) = J_ ^ C , eg 
cos (Aa) E eq C L 
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Figure 14: Experimental results. Ten consecutive figure-eight paths in the (X, Y, Z) frame (solid lines). The dashed lines 
represent the wind window projected on the (X, Y) plane as well as points with constant <fi values and 9 spanning the 
interval [0, 7r/2]. Employed wing: 9 m 2 . Guidance parameters: 6+ = 6- = 0.35 rad, 0_ = — 0.2rad, <fr + = 0.2 rad; 
oj~ = 0.25 Hz. 



where E eq is the equivalent efficiency of the wing. By the equation above we can see that Aa is small for a reasonable wing efficiency 
of 4-6. Measurements of test flights have shown that Aa < 0.3 rad, and most of times Aa < 0.2 rad. Therefore we can linearize (126) 
to get 



Moreover, on the basis of Assumption [2] we can simplify also r\ in {8} as r\ = Aa-i/j. Finally, by Assumption [TJ the component of 
absolute wind perpendicular to the tether is zero, and the velocity angle 7 ~ £. Hence, the aerodynamic force components given by ((6) 
and (7} can be simplified as follows: 




(27) 




(28) 



By using {25} and f2at-<!2bt with {28} we obtain: 




(29) 



Following the same argument by which {27} is derived, we also obtain: 




(30) 



Considering {5}, \v\ is computed as: 




(31) 



which is, by the definition of 7 il lb . also equal to 



\v I = r 9 cos (7) + r cos (8) <f> sin (7) . 



(32) 
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Figure 15: Experimental results. Figure-eight paths flown with different wind speeds. The solid line corresponds to a 
series of consecutive figure-eights flown with average wind speed of 4.5 m/s, the dashed line with average wind speed of 
2.2 m/s. Employed wing: 9 m 2 . Guidance parameters: 9 + = 6- = 0.8 rad, 0_ = — 0.4rad, <fi + — 0.4 rad; u) 1 = 0.25 Hz. 

Finally, by combining d29t with J27t and d30t- d32b . and considering the linearization of (|9) by Assumption[_] we get our result: 

7 = ^(l + 4) 2 |^ + ^|p. +sin ^. 

References 

[1] Ampyx power website, http://www.ampyxpower.com/. 

[2] Windlift website, http://www.windlift.com/ 

[3] Kitenergy website, http://www.kitenergy.net/ 

[4] Enerjite Gmbh website, http://www.enerkite.com/. 

[5] EISG project "Autonomous flexible wings for high-altitude wind energy generation", experimental test movie, August 2012. 
Available on-line: http://lorenzofagiano.altervista.org/movies/EISG_UCSB_auto_wing.mp4 

[6] F. Amato. Quadratic stability. In Robust Control of Linear Systems Subject to Uncertain Time-Varying Parameters, volume 325 
of Lecture Notes in Control and Information Sciences, pages 31-92. Springer Berlin / Heidelberg, 2006. 

[7] I. Argatov, P. Rautakorpi, and R. Silvennoinen. Estimation of the mechanical energy output of the kite wind generator. Renewable 
Energy, 34:1525, 2009. 

[8] J. H. Baayen and W. J. Ockels. Tracking control with adaption of kites. JET Control Theory and Applications, 6(2): 182-191, 
2012. 

[9] M. Canale, L. Fagiano, and M. Milanese. Power kites for wind energy generation. IEEE Control Systems Magazine, 27(6):25-38, 
December 2007. 

[10] M. Canale, L. Fagiano, and M. Milanese. High altitude wind energy generation using controlled power kites. IEEE Transactions 
on Control Systems Technology, 18(2):279 -293, mar. 2010. 

[1 1] Michael Erhard and Hans Strauch. Control of towing kites for seagoing vessels. arXiv, abs/1202.3641, 2012. 

[12] L. Fagiano, K. Huynh, B. Bamieh, and M. Khammash. On sensor fusion for airborne wind energy systems. arXiv, 1211.5060, 
2012. 

[13] L. Fagiano and M. Milanese. Airborne wind energy: an overview. In American Control Conference 2012, pages 3132-3143, 
Montreal, Canada, 2012. 

[14] L. Fagiano, M. Milanese, and D. Piga. High-altitude wind power generation. IEEE Transactions on Energy Conversion, 25(1): 168 
-180, mar. 2010. 

[15] L. Fagiano, M. Milanese, and D. Piga. Optimization of airborne wind energy generators. International Journal of Robust and 
Nonlinear Controll, 22(18):2055-2083, 2011. 




19 



1 

0.9 
0.8 
0.7 

2, 0.6 

a: 

0.5 
0.4 
0.3 
0.2 



Figure 16: Experimental results. Figure-eight paths obtained with different values of the cutoff frequency w 7 . The 
solid line corresponds to u 7 = 0.25 Hz, the dashed one to w 7 = 1Hz. Employed wing: 9 m 2 . Guidance parameters: 
0+ = 6- = 0.55 rad, 0_ = -0.2 rad, 0+ = 0.2 rad. 
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Figure 17: Experimental results. Single figure-eight path obtained during automatic test flights with about 3.1 m/s wind 
speed. From top to bottom: flying path in (<p, 9) coordinates, course of the total force acting on the lines (solid line) and 
of the forces acting on the left (dotted), right (dash-dot) and center (dashed) lines, course of the velocity angle 7 (solid 
line) and reference velocity angle 7 re f (dashed). Employed wing: 12 m 2 . Guidance parameters: 9 + = 9- = 0.4 rad, 
(f>- = -0.2 rad, </>+ = 0.2 rad; w 7 = 0.25 Hz. 
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Figure 18: Experimental results. Single figure-eight path obtained during automatic test flights with about 2.2 m/s wind 
speed. From top to bottom: flying path in (<fi, 9) coordinates, course of the total force acting on the lines (solid line) and 
of the forces acting on the left (dotted), right (dash-dot) and center (dashed) lines, course of the velocity angle 7 (solid 
line) and reference velocity angle 7 re f (dashed). Employed wing: 6 m 2 . Guidance parameters: 9+ = 9- = 0.35 rad, 
(f>- = -0.2 rad, (p+ = 0.2 rad; w 7 = 0.25 Hz. 
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Figure 19: Experimental results. Comparison between the measured values of total line force as a function of the quantity 
E eq ^1 + -]gr\ (cos6 cos{<f))\W\^ (gray dots) and the theoretical linear relationship given by the gain pC ^ — as per 

the results of (T31 (solid lines), (a) Airush One® 6 kite, (b) Airush One® 9 kite, (a) Airush One® 12 kite. The lumped 
parameters for the kites are reported in Table [T] 
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